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We study how the proximity to an integrable point or to localization as one approaches the atomic limit, as 
well as the mixing of symmetries in the chaotic domain, may affect the onset of thermalization in finite one- 
dimensional systems. We consider systems of hard-core bosons at half-filling with nearest neighbor hopping 
and interaction, and next-nearest neighbor interaction. The latter breaks integrability and induces a ground- 
state superfluid to insulator transition. By full exact diagonalization, we study chaos indicators and few-body 
observables. We show that when different symmetry sectors are mixed, chaos indicators associated with the 
eigenvectors, contrary to those related to the eigenvalues, capture the onset of chaos. The results for the com- 
plexity of the eigenvectors and for the expectation values of few-body observables confirm the validity of the 
eigenstate thermalization hypothesis in the chaotic regime, and therefore the occurrence of thermalization. We 
also study the properties of the off-diagonal matrix elements of few-body observables in relation to the transition 
from integrability to chaos and from chaos to localization. 
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I. INTRODUCTION 

The usual approach to the study of systems with complex 
energy spectra, such as nuclei, atoms, molecules, and quan- 
tum dots, is via random matrices. These are matrices filled 
with random numbers whose sole restriction is to satisfy the 
symmetries of the system under investigation 111-0]. Time- 
reversal invariant systems with rotational symmetry, for in- 
stance, are described by the so-called Gaussian orthogonal 
ensembles (GOEs), which consist of ensembles of real sym- 
metric random matrices. The range of applicability of random 
matrix theory (RMT) was further extended after the connec- 
tion with classical chaos became established. It was verified 
that the spectra of quantum systems that behave chaotically 
in the classical limit show the same fluctuation properties ob- 
tained with ensembles of random matrices. This observation 
was stated in the form of a conjecture ^ and initiated the field 
of quantum chaos. 

The most commonly used quantities to identify the onset of 
quantum chaos are based on the eigenvalues of the Hamilto- 
nian describing the system under investigation, however the 
structures of the eigenvectors have also a very important role 
ll6l- fl0ll . The eigenvectors of a system whose classical coun- 
terpart is chaotic are expected to be maximally delocalized. 
According to Berry's conjecture ITTl [l2ll . the eigenfunctions 
become superpositions of plane waves with random phases 
and Gaussian random amplitudes. 

Berry's conjecture has been connected to the problem of 
thermalization in isolated quantum systems lisi [13U16I1 . As 
stated in Ref. 13, "a bounded, isolated quantum system of 
many particles in a specific initial state will approach thermal 
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equilibrium if the energy eigenfunctions which are superposed 
to form that state obey Berry's conjecture". It is possible to 
show that such eigenstates lead to the appropriate (Maxwell- 
Boltzmann, Bose-Einstein, or Fermi-Dirac) distribution for 
the momentum of the particles in the system ifTill . In this sce- 
nario, eigenstate expectation values (EEVs) do not fluctuate 
between eigenstates that are close in energy and hence they 
coincide with the microcanonical average. The latter became 
known as the eigenstate thermalization hypothesis (ETH). 

The interest in the problem of thermalization, and in the 
dynamics of isolated quantum systems far from equilibrium 
in general, was recently boosted by experiments with ultra- 
cold atoms in optical lattices. In particular, the antagonistic 
results obtained with a bosonic gas in quasi-one-dimensional 
geometries, where thermalization was inferred to occur in one 
experiment 1 1 7| but was not observed in another one ifisll . 
motivated various theoretical studies of nonintegrable one- 
dimensional (ID) quantum systems after a quench. A spe- 
cial property of the ID systems that have been analyzed is the 
possibility to reach both integrable and nonintegrable regimes 
by adjusting parameters of the Hamiltonian. It was verified 
that close to the integrable point, ETH ceases to be valid and 
thermalization does not happen lfT9ll20ll . But the absence of 
thermalization has also been linked to other factors, such as 
the effects of particle statistics in finite systems 1I20I1 and the 
proximity of the energy of the initial state to the energy of the 
ground state of the system after the quench lfT9l - l2lll . 

A close inspection of the static properties of the models be- 
ing studied can anticipate the results for the dynamics ll22ll . In 
finite ID lattices, the integrable-chaos transition for fermions 
has been shown to require larger integrability-breaking terms 
than for bosons, which can explain the lack of thermalization 
of the former in certain regimes ll22ll . Models describing re- 
alistic systems involve only few-body interactions, therefore 
random matrices need to be substituted by banded matrices 
lfToll23l - l26ll . In clean systems, such as the ones involved in re- 
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cent studies lfT9[|20ll22ll27ll . these sparse matrices do not even 
contain random elements. Full random matrices and banded 
matrices may show similar spectral statistics, but they differ in 
terms of level density, the first showing a semicircular spec- 
trum and the latter a Gaussian spectrum, and in terms of eigen- 
states. Contrary to full random matrices, where all eigenstates 
are random vectors, in the case of finite range interactions, 
chaos develops only away from the edges of the spectrum, 
so it is only there that the eigenstates can satisfy ETH. This 
explains why nonequilibrium initial states with energy close 
to the borders of the spectrum are not expected to thermal- 
ize B . 

Further factors that have been associated with the absence 
of thermalization in finite ID systems are the opening of a gap 
as one crosses a superfluid to insulator transition |!28:1|, and the 
existence of "rare" states, which for nonintegrable systems 
have been speculated to persist in the thermodynamic limit 
ll29ll . In general, the question of thermalization as one crosses 
a superflu id ( metal) to insulator transition has attracted a lot of 
attention ll2ll l27l - l3lll . We have recently argued that thermal- 
ization does happen in the gapped side of the phase diagram, 
and that as one increases the system size it occurs deeper into 
that side IztIi . We did not find evidence of the existence of 
rare states in those systems ll27ll . Thermalization ceased to 
occur only when the system approached the atomic limit and 
the eigenstates started to localize in the momentum basis. 

In the present work, we further analyze the issue of thermal- 
ization in systems that approach a localization regime close to 
the atomic limit. As in Ref. |27], localization refers here to 
the broad notion of contraction of the eigenstates in a partic- 
ular basis set, instead of the more specific concept of spatial 
localization due to disorder. Disorder is absent in the systems 
that we consider. In comparison to Ref. IztIi . an extra com- 
plication is added to our studies: when analyzing the observ- 
ables of interest, some discrete symmetries are not removed. 
We then address the role of such symmetries in their static 
and dynamical properties. It is well known that the mixing of 
symmetries may conceal key features of the chaotic regime, 
such as level repulsion ll32ll33ll . Could it affect also the va- 
lidity of the ETH? We show that the structure of eigenvectors 
that belong to different subspaces remain very similar in the 
chaotic region. As a result, EEVs do not fluctuate and the 
ETH continues to be valid. 

We focus on ID systems of hard-core bosons (HCBs) at 
half-filling, with nearest neighbor (NN) hopping (t) and in- 
teraction (V) and with next-nearest neighbor (NNN) interac- 
tion (V). In the absence of NNN interaction, the model is 
integrable; we study how integrability is broken by V'. In 
addition, when V' ^ t,V a transition to localization in the 
momentum basis starts to take place and is reflected in the in- 
verse participation ratio (IPR) of the eigenstates of the Hamil- 
tonian. We study how this second transition affects the EEVs 
of various few -body observables and their off-diagonal matrix 
elements. EEVs and the off-diagonal elements of the observ- 
ables are related to the occurrence of thermalization and to the 
time evolution of the system after a quench, respectively. 

The paper is organized as follows. Sectionlllldescribes the 
model Hamiltonian studied and its symmetries. SectionHIJan- 



alyzes the integrable-chaos transition based on various chaos 
indicators. The eigenstate expectation values of different ob- 
servables and the comparison with the microcanonical aver- 
ages are shown in Sec. |IVl Section [V] is devoted to study- 
ing the behavior of the off-diagonal elements of few-body ob- 
servables in the eigenstates of the Hamiltonian. Concluding 
remarks are presented in Sec. |Vll Further illustrations about 
delocalization measures and observables are provided in the 
Appendix. 

II. SYSTEM MODEL 

As mentioned in the introduction, we study a ID HCB 
model with NN hopping t and interaction V, and NNN in- 
teraction V'. The Hamiltonian is given by 

L 

Ht = Y,{-t{b}M+i+ii-c) (1) 

i=l 

+V (n, - i) - + r (n, - (n.^, - } . 

where L is the size of the chain, bi (hi) is the bosonic annihi- 
lation (creation) operator on site i and fii = is the boson 
local density operator Hard-core bosons are not allowed to 
occupy the same site, so &| = = 0. 

Hamiltonian ([TJ conserves the total number of particles 
Ni, and is translational invariant. It consists of independent 
blocks, where each one is associated with a value of Ni, and a 
total momentum k. Here we study chains with an even num- 
ber of sites and at half-filling, Nf, = L/2, and consider all 
values of k, from to i/2. At half-filling, other symmetries 
are found; particle-hole exists for all fc's and parity appears 
only for k = 0, L/2. We perform full exact diagonalization 
of each fc-sector separately for chains of 18, 20 and 22 sites. 
The dimension Dk of each fc-sector is given in Table|T] Notice 
that we need to take into account also the double multiplicity 
of the eigenstates belonging tofc=:l,...L/2 — 1. The largest 
total Hilbert space considered has dimension D = 705 432. 



TABLE I: Dimension of fc-sectors 



L = 18 
Dk 


= 0,9 
2704 


k = 1,2,4,5,7,8 
2700 


fc = 3,6 
2703 




L = 20 
Dk 


fc = 0, 10 
9252 


= 1,3,7,9 
9225 


fc = 2,4,6,8 
9250 


fc = 5 
9226 


L = 22 
Dk 


fc = 0, 11 
32066 


other fc's 
32065 







In what follows, t ^ 1 (h = I) sets the energy scale and 
the interactions are repulsive, V, V > 0. We fix = 1 and 
vary V' from to 10. The system is integrable when V' = 0, 
while the addition of NNN interaction may lead to the onset 
of chaos. Moreover, there is a critical value of the NNN inter- 
action, V/ = 2, below which the ground state is a gapless su- 
perfluid and above which it becomes a gapped insulator Isill . 
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A small bond-ordered phase develops around = 2 Issll . 
When V' ^ t, V, the system approaches the atomic limit. 
Due to the translational invariance of model ([TJ, the eigen- 
states of the Hamiltonian approach the eigenstates of the total 
momentum operator. 

III. QUANTUM CHAOS INDICATORS 

Notions of phase-space trajectory and Lyapunov exponent, 
which are used to distinguish regular from chaotic motion 
in classical mechanics, have no meaning in the quantum do- 
main. Nevertheless, criteria exist to separate quantum systems 
whose classical counterpart are chaotic from those whose 
classical counterpart are regular. Signatures of quantum chaos 
are obtained from the eigenvalues and from the eigenvectors 
of the Hamiltonian. 



A. Spectral observables 

Spectral observables, such as level spacing distribution, 
level number variance, and spectral rigidity are intrinsic in- 
dicators of the integrable-chaos transition However, a 
main disadvantage associated with the computation of these 
quantities is the need to identify and separate all symmetry 
sectors of the system. It is only after the separation that the 
spectrum may be unfolded and the analysis carried out. 

The distribution of spacings s of neighboring energy levels 
is the most frequently used observable to study short-range 
fluctuations in the spectrum Quantum levels of inte- 

grable systems are not prohibited from crossing and the dis- 
tribution is Poissonian, 

Pp{s) = cxp(-s). (2) 

In non-integrable systems, crossings are avoided and the level 
spacing distribution is given by the Wigner-Dyson distribu- 
tion, as predicted by RMT. The form of the Wigner-Dyson 
distribution depends on the symmetry properties of the Hamil- 
tonian. Ensembles of random matrices with time reversal in- 
variance and rotational symmetry, the GOEs, lead to 

TTS f IT \ 

Pw d{s) ^ —cyi^i — ^j- (3) 

The same distribution form is expected for Hamiltonian ([T]i in 
the chaotic regime, even though Hb has only two-body inter- 
actions and does not contain random elements. 

The top panels of Fig. [T] depict the level spacing distribu- 
tions for Hamiltonian ^ when L = 22. P{s) is computed 
for each fc-sector separately and the results are then averaged 
between fc = and k = L/2, and between the rest of the 
/c-sectors. The decision to perform two different averages is 
made because inside each /c-sector, particle-hole symmetry is 
present for all fc's, but parity exists only for k = 0, L/2; in 
addition, all the sectors in each average behave very similarly. 
The top panels in Fig. [T] clearly show that the distribution of 
spacings never becomes equal to Pwd{s), instead they show 



an intermediate behavior between Pp{s) and Pwd{s), even 
though we do expect these systems to be chaotic away from 
V' = 0. This occurs because we have not separated the sub- 
spaces according to all symmetries. The mixing of the re- 
maining symmetries in each fc-sector obscures the effects of 
level repulsion lf36ll . Notice that, as expected from the amount 
of remaining symmetries, sectors fc = and fc ~ i/2 are the 
ones further away from a Pwd{s) distribution. The overall 
behavior of P{s) in these systems is certainly in contrast with 
the results presented in Refs. ll22ll27ll for other ID chains, and 
shows that the presence of as many as one discrete symmetry 
may hinder the signatures of quantum chaos. 
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FIG. 1 : (Color online.) Top panels: Level spacing distributions. Re- 
sults are shown for the average between sectors k = Q and k = L/2 
[black solid line], the average between all other fc-sectors [light (red) 
solid line], the Poisson distribution, and the Wigner-Dyson (W-D) 
distribution. Bottom panels: Inverse participation ratio in momen- 
tum space for sectors fc = and fc = L/2 [black dots] and all other 
fc-sectors [light (red) dots]. In all cases the system size is L = 22. 



B. Delocalization measures 

Quantities that focus on the eigenvectors, as delocalization 
measures are not intrinsic indicators of the integrable- 

chaos transition, since they depend on the basis in which the 
computations are performed. But, contrary to spectral observ- 
ables, we show here that these quantities do not necessarily 
require a separated analysis of each symmetry sector. 

The degree of delocalization of individual eigenvectors may 
be measured, for example, with the inverse participation ra- 
tio (IPR), denoted here by ^, or the information (Shannon) 
entropy S . The former is also sometimes referred to 
as number of principal components (NFC). For an eigenstate 
\'^a) of Hamiltonian ([T]i written in the basis vectors as 
Y^a) ~ Sjii ^il^j), IPR and S are respectively given by 
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and 

S^^-J2\c>j'\n\c'j'. (5) 

The above quantities measure the number of basis vectors that 
contribute to each eigenstate, that is, how spread each state is 
in the chosen basis. 

The choice of basis is usually determined by the infor- 
mation one is after and by possible computational limita- 
tions. Here, we consider the momentum basis, given that for 
large values of V' the system exhibits localization in /c-space. 
This is an interesting effect that results from approaching the 
atomic limit while imposing translational symmetry. Other 
relevant bases include the mean-field basis [8], which corre- 
sponds to choosing the eigenstates of the integrable Hamilto- 
nian {V' = 0) as a basis and therefore captures localization as 
y — > ll22i I27I1 . and the site basis, which is meaningful in 
studies of spatial localization. 

GOEs lead to extreme delocalization, their eigenvectors 
are random vectors where the amplitudes c" are indepen- 
dent random numbers. The average over the ensemble gives 
S'^"^ ~ ln(0.48i:»fc) and _ ^^^^^^ Hamilto- 

nian ^ has only two-body interactions, the eigenstates of our 
system in the chaotic limit may approach the GOE result only 
away from the edges of the spectrum ll24l |26[ Issll . 

The bottom panels of Fig. [T] show IPR in the fc-basis (S^k) 
for all fc-sectors for the same values of V' used in the study of 
the level spacing distributions (top panels). While the results 
for P{s) hardly change with V', three different regimes can 
be singled out from the behavior of ^fc. (i) When V' —5- 0, 
the values of ^fc fluctuate considerably for states very close 
in energy, which agrees with our expectations for a system in 
the integrable regime ll22ll . (ii) For intermediate values of V' 
[0 < V^' < 5 for L = 22], becomes a smooth function of 
energy, which suggests the crossover to chaos ll22ll . (iii) At 
large values of V', energy bands are created and the once 
again fluctuates considerably. We also notice that in the sce- 
nario (iii), £,k decreases significantly, signaling localization in 
/c-space. The two transitions IztIi . from integrability to chaos 
as V' increases from zero and from chaos to localization in 
/c-space as V' — > cxd, are therefore clearly captured by 
despite the inclusion of eigenstates from different subspaces. 

Two separated curves are clearly distinguished in the bot- 
tom panels of Fig.[T]when V' ~ 1, 2, 3, 4. This is caused by 
two combined factors. First, particle-hole symmetry exists for 
all fc-sectors, but parity is only present for fc = and k = L/2. 
Thus, the eigenstates from the two latter subspaces cannot 
spread as much as those pertaining to fc = 1,2, ...L/2 — 1 
and so exhibit smaller values of £^k- Second, the structures of 
the eigenstates from different fc-sectors, containing the same 
number of internal symmetries, are very similar and do not 
fluctuate in the chaotic region. As a result, the two domains, 
regular with large fluctuations and chaotic with no fluctua- 
tions, are well distinguished. This explains why the analysis 
of is an efficient way to detect the transition to chaos even 
when we do not separate the eigenstates according to all of 
their symmetry sectors. 



With increasing system size, we find that chaotic behavior 
beyond V/ can be observed for larger values of V'. This is 
shown in the Appendix (A.l), where we argue that similar 
conclusions, as the ones presented here with (^k, are reached 
with the Shannon entropy. 

C. Structural entropy 

Further information about the structure of the eigenvectors 
may be obtained with the so-called structural entropy, which 
is defined as ll39l|40|l 

S,tr = S-lne (6) 

Sstr contains the contribution to the information entropy which 
is not found in the IPR. It is an attempt to better distinguish 
states that may have similar levels of delocalization, but dif- 
ferent structures. In the case of a GOE, the states are uniform 
and SH?^ « 0.3646. 
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FIG. 2: (Color online.) Structural entropy in the momentum basis vs 
energy per site for all fc-sectors. Black dots: L = 22, light (green) 
dots: L — 20. The values of V' are indicated in the panels. The solid 
line correspond to f» 0.3646. 

Since Ssu does not aim at measuring the actual extension of 
the eigenstates, but instead at capturing their structures, vec- 
tors belonging to different symmetry sectors may be analyzed 
on a par with each other, even when they have different levels 
of delocalization. In Fig.|2]we show Sstr for all fc-sectors and 
for two system sizes, L ~ 20 and L = 22. The results for all 
sectors are remarkably similar and superpose each other The 
plots reveal again the same three regimes identified in Fig. [T] 
(i) For V' -> 0, eigenstates very close in energy have different 
levels of complexity, as typical of integrable systems, (ii) The 
structures of the states become comparable to random vectors 
in the middle of the spectrum when chaos is reached [when 
< V' < 5 for L = 22]. The two-body interactions are 
responsible for the bowl-shaped curve and the fluctuations at 
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the edges of the spectrum, (iii) As the system moves to local- 
ization in the fc-basis, for large V', energy bands accompanied 
by large fluctuations of the values of Ssu become evident. The 
fact that the analysis of Ssu does not require the identifica- 
tion and separation of symmetry sectors supports our claim 
that quantities associated with the eigenvectors may, in many 
instances, be better suited than spectral observables for study- 
ing the integrable-chaos transition, especially when unknown 
symmetries may be present. 

In Fig.m we present results for two system sizes, L = 20 
and L = 22. The results are very similar, but Ssu- for L = 20 
exhibits larger fluctuations, particularly in the chaotic region, 
i.e. fluctuations decrease in the chaotic regime as L increases. 
Also, for large values of V', where eigenstates are grouped in 
bands with similar energies, Ssn- shows that those bands shift 
as the system size increases. 
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FIG. 3: (Color online.) Standard deviation of Sstr vs V' for eigen- 
states in the middle of the spectrum with energies varying from -5 to 
5. The inset shows the gap (between the lowest energy states) times 
L vs V' . Results are shown for lattices with L — 18, L = 20, and 
L = 22 sites. 

In Fig.[3]we compare the standard deviation of the structural 
entropy, 

a(S.„) = y^(S?J-(S.tr)2, (7) 

for states in the middle of the spectrum for different lattice 
sizes. As V' increases from zero, <t captures the two tran- 
sitions of model This is particularly visible for L = 
22, where large fluctuations appear in the integrable domain 
{V — !> 0) and in the localization regime {V' oo), while 
small fluctuations are associated with the onset of chaos. 

In relation to the low energy behavior of these systems, the 
opening of the gap A between the ground state and the first 
excited state, signaling the onset of the superfluid to insulator 
transition (the ground state becomes four-fold degenerate in 
the insulating side) is well illustrated by the curve for i = 20 
in the inset of Fig. |3] [boundary effects conceal the transition 
for L = 18 and 22]. By comparing the value of V' for the 
chaos-localization transition with the value for the superfluid- 
insulator transition, it becomes clear that an overlap between 
chaotic regime and gapped phase exists for the finite systems 
considered here. In the case of L = 20, for instance, the 



opening of the gap is already evident when V ^ 3 (see inset 
of Fig. [3), while the formation of energy bands followed by 
the localization in the momentum basis requires V > A (see 
Figs.|2l[3]and[l2li. 

The dispersion in the main panel of Fig. [3] makes evident 
also the dependence of the results on L. Larger systems imply 
smaller fluctuations. Moreover, as L increases, smaller val- 
ues of V' already lead to the first transition from integrability 
to chaos and larger values of V are required for the second 
transition from chaos to localization. The shift of V' to larger 
values for the second transition shows that, as the system in- 
creases, chaoticity appears deeper into the gapped phase. In 
the thermodynamic limit, one may even speculate that any 
V Q might suffice to guarantee the chaoticity of the system. 
These results reinforce the claim that the superfluid-insulator 
transition does not affect the behavior of the bulk of the eigen- 
state of the Hamiltonian ll27ll . 

The uniformization of the eigenvectors in the chaotic 
regime has been manifested in our studies of IPR, Ssa and S 
in Figs. [ni2l 1311 II and[T2l The results prompt us to advocate, 
as in Ref. 11411 and references therein, that in certain situations 
quantities to measure the complexity of the eigenstates may be 
better indicators of quantum chaos than spectral observables. 
The analysis of the eigenstates hints also on what to expect in 
terms of thermalization. Thermalization has since long been 
associated with chaos and ergodicity. At the classical level, 
the idea is well established [|4111 . while in the quantum domain 
the connection is based on a hypothesis, the ETH. Accord- 
ing to the ETH 13, the eigenstate expectation values (EEVs) 
of few-body observables do not fluctuate between eigenstates 
that are close in energy and hence they coincide with the mi- 
crocanonical average. This reflects the fact that in the chaotic 
regime the structure of the eigenstates in a small interval of 
energy may be thought as equivalent. The smooth behavior of 
EEVs with energy, which is achieved in the chaotic domain, 
is discussed and illustrated in the next section. 

IV. FEW-BODY OBSERVABLES 

Here, we provide numerical support for the connection be- 
tween the ETH and quantum chaos. This is done based on the 
analysis of the EEVs of four different observables: 

(i) the kinetic energy, 

k = -tY,{bYh+i+W.c}j, (8) 

i 

(ii) the interaction energy, 

(iii) the momentum distribution function, 

i,3 
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(iv) and the density-density correlation structure factor. 



'niUj. 



(11) 



Since the operator for the total number of bosons commutes 
with the Hamiltonian, the expectation value of N{k = 0) is 
simply {N{k = 0)) ~ ^b/-^- This value is set to zero in 
what follows. K and n{k) are one-body observables, local 
and non-local, respectively; while / and N{k) are two-body 
observables, local and non-local, respectively. K and n{k) are 
routinely measured in cold gases experiments. 
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FIG. 4: (Color online.) Eigenstate expectation values of K and / vs 
energy per site for the full spectrum, including all momentum sec- 
tors. Black dots; L = 22, light (green) dots: L = 20. 

Figures |4] and |5] show the EEVs for the four observables 
defined above. The results parallel the findings for the eigen- 
states obtained in the previous section. As V' increases and 
one departs from integrability, the fluctuations are signifi- 
cantly reduced away from the borders of the spectrum, and 
ETH becomes valid. It is remarkable that despite the inclu- 
sion of EEVs for all disconnected fc-sectors, the results are still 
very similar for eigenstates that are close in energy lfl9il20ll . 
Contrary to the eigenstates, which showed at least a differ- 
ence in the level of delocalization depending on their fc-sector, 
states for k = being less spread than for the other fc's 

(cf. bottom of Fig. [T]!, the behavior of the EEVs in all fc- 
sectors is very similar Hence, it is not surprising that if some 
discrete symmetries are not accounted for when diagonalizing 
the Hamiltonian, ETH will still be valid in the chaotic regime. 
No separation of the EEVs occurs for different symmetry sec- 
tors. 

The smooth behavior of EEVs with energy, which is char- 
acteristic of the chaotic domain, continues to hold beyond the 
superfluid-insulator transition (compare Figs. |4] and |5] with the 
inset of Fig.O, further confirming that the latter is irrelevant 
for the discussion of the validity of ETH. By increasing V' 
even further, the eigenstates finally begin to localize in fc- 



space and large variations of the EEVs for states close in en- 
ergy reappear. In this limit, the separation of the expectation 
values into energy bands becomes evident. 
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FIG. 5: (Color online.) Eigenstate expectation values of n{k) and 
N{k) vs energy per site for the full spectrum, including all momen- 
tum sectors. Black dots: L — 22, light (green) dots: L — 20. 

The comparison between the results for different system 
sizes in Figs |4] and |5] L = 20 and L ~ 22, shows that (i) 
as the system size increases, the fluctuations between EEVs 
of states that are close in energy decrease in the chaotic re- 
gion, and (ii) in the localization regime, which is reached for 
large V', the position of the energy bands for L = 20 and 
L ~ 22 do not coincide. The comparison also reinforces the 
disconnection between the behavior of low energy states and 
the bulk of states. As seen in the inset of Fig. [3] a gap opened 
for L ~ 20, but boundary effects prevented it in L = 22. This 
difference has no consequences in the results for EEVs, which 
are very similar for both L's. 



A. Eigenstate thermalization hypothesis 

Strong evidence of the validity of the ETH is established 
once EEVs are seen to be very similar between eigenstates 
that are close in energy. This, in turn, implies that thermal 
averages and the EEVs will be also very similar The analy- 
sis above shows that this should occur in the chaotic regime. 
To quantify this statement for finite systems, we compute the 
deviation of the EEVs for an observable O with respect to the 
microcanonical result (A""^), defined as 



(12) 



In Eq. ( IT2] |. the sum runs over the microcanonical window, 
Oaa are the EEVs of the operator O, and the microcanonical 
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expectation values Omic are obtained from 
1 



^mic — 



E.AE 



E 



where Ne,ae is the number of energy eigenstates with energy 
in the window [E — AE, E + AE]. In what follows, we will 
also refer to A™"^© as the average fluctuations of the EEVs. 




FIG. 6: (Color online.) Average relative deviation of eigenstate 
expectation values with respect to the microcanonical result as a 
function of V' . The average is performed over all eigenstates (in- 
cluding all momentum sectors) that lie within the window [E — 
AE,E + AE], with AE = 0.4. Results are shown for lattices 
with L — 18, L — 20, and L — 22 sites. The effective tempera- 
ture is T = 5, which for L = 22 corresponds to E — —2.9680 for 
V = 0, E = -4.3126 for V ^ 2, E = -7.7462 for V = 4, 
E = -13.1636 for V = G, E = -20.4001 for V = 8, and 
E = -29.2205 for V = 10. 

Figure |6] shows the relative deviations A™"^/^, A™"^/, 
A™'^n(fc = 0), and A™"^ (A; = tt) averaged over all momen- 
tum sectors and for all eigenstates that lie within a window 
[E - AE,E + AE], where AE = 0.4. While the results 
should not depend on the exact value of AE around a reason- 
able choice, the selection becomes subtle for large values of 
V', where energy bands are formed and the number of states 
for small energy windows decay significantly. Our choice was 
made to guarantee that all fc-sectors, for all system sizes and 
for all values of V', have a sufficiently large number of eigen- 
states in [E — AE,E + AE]. [For more discussion of the 
dependence of the results on AE, see the Appendix (A. 2)]. 
The value of E is selected according to the effective temper- 
ature T that we chose to study. Performing the analysis in 
terms of a single temperature allows for a fair comparison of 
all systems sizes and values of V'. The effective temperature, 
Ta of an eigenstate with energy Ea is defined as 



Ea — 



where 



Z = Tr{e-^/^°} 



Above, H is Hamiltonian Z is the partition function with 
Boltzmann constant ~ I, and the trace is performed over 
the full spectrum. 

As seen in Fig.|6] the average fluctuations of the EEVs for 
all observables considered decrease as V' increases from zero 
and the integrable-chaos transition takes place, which goes 
along with the validity of the ETH in the chaotic domain. The 
dependence of the results on system size is also clear; the 
average fluctuations decrease for larger systems. In addition, 
the width of the interval of values of V' for which the EEVs 
approach the thermal averages can, in general, be seen to in- 
crease with L, which brings the validity of ETH deeper into 
the gapped phase. On the other hand, beyond the chaotic do- 
main, as y — > oo and the system approaches the atomic limit, 
large fluctuations reappear As the system starts to localize 
in the momentum basis and energy bands are formed (when 
V^' > 6 for L 22), A'^'K and A™'=n(fc = 0) become 
even larger than in the integrable regime [cf. panels (a) and 
(c)]. This is a consequence of the large fluctuations that occur 
especially for the kinetic energy and the momentum distribu- 
tion function (see Figs.|4]and|5]i in the windows of eigenstates 
associated with the chosen effective temperature T = 5. 

As discussed in Ref. ll22ll . the proximity to the ground 
state prevents thermalization in systems with few-body in- 
teractions, even in the chaotic domain, since chaos does not 
develop at the edges of the spectrum. For the systems con- 
sidered here, the presence of a gap is an additional hindering 
factor for the thermalization of nonequilibrium initial states 
with low energies. This is because far from the ground state 
now means that the energy of the time-evolving state has to 
be greater than the energy of the first excited state, which is 
determined by the gap. This implies that the minimal effec- 
tive temperature at which thermalization will occur increases 
as the gap in the system increases. 

The study of the average fluctuations of the EEVs as a func- 
tion of temperature further supports the conclusions above. 
In Fig. |2l we present results for A™'^n(fc = 0) vs T for ten 
different values of V' and for temperatures T < 10. As 
V^' — > and we approach the integrable regime, large val- 
ues of A™"^n(A; = 0) appear for all temperatures considered, 
[cf. panel (a)]. Contrary to that, in the chaotic domain, large 
values of A'™'^n(fc = 0) are restricted to low temperatures, 
while at large T, A'"'^n{k ~ 0) saturates at small values [cf. 
panels (b) and (c)]. This corroborates our statements that the 
validity of ETH goes hand in hand with the onset of chaos and 
holds away from the edges of the spectrum. Far from chaotic- 
ity, when V' oo, large fluctuations are seen for various 
temperatures, and the peaks of A™"^n(fc = 0), associated with 
the energy bands, move in temperature as L increases. 

It has been discussed in Ref. ll29ll that, for local observ- 
ables, the deviation of the EEVs from the microcanonical 
average [given by Eq. ( [T2b 1 vanishes as the system size in- 
creases. This result is independent of whether the system 
is integrable or not. In our figures in this section, we have 
clearly shown that, for any given system size, the deviation 
of the EEVs from the microcanonical average is, in general, 
larger away from the chaotic regime, no matter whether the 
observable is local or nonlocal. How those fluctuations vanish 
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FIG. 7: (Color online.) Average relative deviation of eigenstate 
expectation values of h{k) with respect to the microcanonical pre- 
diction as a function of the effective temperature T of the eigen- 
states. Panels (a) and (d) depict results for L = 20 and L = 22. 
In panels (b) and (c) only L = 22 is shown, since the results for 
L = 20 are very similar. The average is performed over all eigen- 
states (including all momentum sectors) that lie within the window 
[E - AE, E + AE], with AE = 0.1. 
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FIG. 8: (Color online.) Normalized extremal fluctuations of eigen- 
state expectation values as a function of V' . All eigenstates from all 
momentum sectors that lie within the window [E — AE, E + AE], 
with AE = 0.4 are taken into account. Results are shown for lattices 
with L = 18, L = 20, and L = 22 sites. The effective temperature 
is T = 5, which for L = 22 con-esponds to E ^ -2.9680 for 
V ^ 0, E = -4.3126 for V ^ 2, E = -7.7462 for V = 4, 
E = -13.1636 for V ^ 6, E = -20.4001 for V = 8, and 
E = -29.2205 for V' = 10. 



as the system size increases can depend on whether the sys- 
tem is integrable or not, and is something that deserves further 
investigation. In Fig. |6l the deviations of the EEVs from the 
microcanonical result, in particular for the nonlocal observ- 
ables n(k) and N{k), are seen to decrease faster with system 
size in the chaotic regime. 

We should stress, however, that our calculations in the 
chaotic regime not only show that the average deviations of 
EEVs for all our observables decreases as one increases the 
system size, but also that the same occurs with the extremal 
fluctuations of the individual EEVs. This provides a more rig- 
orous test of the validity of the ETH. 

We have studied the normalized extremal fluctuation of an 
observable O, defined as. 



Aro EE 



:0- 



(13) 



The maximum and minimum values of Oaa, maxO and 
minO, are extracted from the same energy window [E — 
AE, E + AE] used to obtain the microcanonical expectation 
value. 

Figure [8] shows the normalized extremal fluctuations 
Af'^K, Af'^I, A™'=n(fc ^ 0), and A™'=iV(fc = tt) for all 
eigenstates from all momentum sectors that lie within a win- 
dow [E - AE, E + AE], where AE = 0.4. This corresponds 
to the worst scenario, where maximum and minimum values 
of Oaa may belong to different fc-sectors. Figure |8] mirrors 
some of the features already seen in Fig. |6] for any given 
system size, as V' increases, the extremal fluctuations for all 
considered observables first decrease as the integrable-chaos 
transition takes place and then increase again as the systems 



approaches the atomic limit. However, in contrast to Fig. |6] 
Fig. [8] makes it evident that the extremal fluctuations of the 
EEVs for all observables decrease with increasing system size 
only in the chaotic region, as expected for the validity of the 
ETH. 



V. PREDICTIONS FOR THE DYNAMICS 

In this section, we discuss what to expect for the time evo- 
lution of an arbitrary initial state under the unitary dynamics 
dictated by Hamiltonian ([T]). 

For an isolated quantum system, the time evolution of an 
initial state \ipini) is determined by 

a 

where \'^a) are the eigenstates of the Hamiltonian and Ca = 
(^alV'mi)- The expectation value of an observable O at time 
t is given by 

{d{t)) ^ {m\o\m) = 5]c:c;3e^(^-^^)*o„^, (u) 

afl 



where 



O. 



Q/3 



{^o.\0\^p) 



are the matrix elements of O in the energy eigenstate basis. 
The infinite time average of the observable corresponds to 



(15) 
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where "diag" stands for diagonal ensemble, that is an ensem- 
ble where each state has weight |C^p ifTel [l9i |20|l . 

Our results so far indicate that Odiag will coincide with Omic, 
that is thermalization will occur in the chaotic regime (where 
ETH is valid), whenever the initial state has an expectation 
value of the energy that is not close to the edges of the spec- 
trum and for a distribution of |Cq | that is sufficiently narrow. 
The latter has been argued to be the case for generic quenches 
lfl6ll . When ETH is not valid, the outcome of Odiaa will de- 
pend on the details of the weights | p and will not be, in 
general, in agreement with the predictions of standard ensem- 
bles of statistical mechanics. 

The question we address here is what to expect for the time 
that will take for the initial state to relax to the diagonal en- 
semble predictions (relaxation time) and for the time fluctua- 
tions that will occur about such an infinite time average. We 
might expect that longer relaxation times, as well as enhanced 
fluctuations after relaxation, should take place close to the in- 
tegrable and localization regimes, and to edges of the spec- 
trum in the chaotic region. These three scenarios may reduce 
the number of states j^'^) with a relevant role in the evolution 
of the initial state and therefore reduce the effects of dephas- 
ing in Eq. (fT4b . Interestingly, in previous works, the relax- 
ation time at and close to integrability has not been found 
to be much different from the one away from integrability 
lfl9[ I20I I27I1 . On the other hand, the approach to localization 
in Ref. iZTIl was shown to substantially increase the relaxation 
time and the time fluctuations after relaxation (see Fig. 3(d) in 
iIztII ). Other factors that may play a role are analyzed in what 
follows. 

A quick relaxation to Odiag and reduced time fluctuations 
require a nondegenerate and incommensurate spectrum, as ex- 
pected for nonintegrable systems. However, as it has been dis- 
cussed in this work, mixing of symmetries may occur even in 
the chaotic domain. In this case, states very close in energy 
may appear and one may wonder if they could slow down the 
dephasing process in Eq. (fl4] i. 

In studies of the unitary dynamics, the system is usually 
taken out of equilibrium by means of a quench. One starts 
with an initial state of a certain Hamiltonian Hini and then 
instantaneo usly changes it to Hfin at time t = 0. In previous 
works lfl9l I20I I27I1 . Hini and Hfin involved the same sym- 
metries, and the initial state was taken from the fc = sec- 
tor, which had an internal remaining symmetry, parity (these 
systems were at 1/3 filling). Surprisingly, the relaxation dy- 
namics in the integrable and near integrable regimes, as well 
as in the chaotic regime, were very similar il9ll20ll27ll . even 
though the level spacing distributions for both domains are 
clearly different. Based on those results, we expect a similar 
behavior for the cases considered in this work, away from the 
localized regime, even if some discrete symmetries remain in 
the fc-sector where the dynamics is performed ll33ll . This is, 
however, another subject that deserves further investigation. 

From Eqs. ( fT4] i and ( flSl l. one realizes that the time fluc- 
tuations of a particular observable O after relaxation can be 



quantified by the expression 

(0(t)) - Jdit)) = E C:C;3e^(^°-^^)*0„^, (16) 

which means that the off-diagonal matrix elements of the ob- 
servable under consideration play a very important role. 
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FIG. 9: (Color online.) Matrix elements of n{k — 0) [top panels] 
and N{k = tt) [bottom panels]; L = 22, k = 0. We select as 
the central state the one that has the closest energy to the energy 
in the canonical ensemble corresponding to an effective temperature 
T — 5.0. For the cases depicted in the figure this corresponds to i5 = 
-2.9680 for V = 0, £ = -3.3748 for V = 1, E ^ -4.3126 
for V ^ 2, E = -13.1636 for V = 6, and E = -29.2205 
for V' = 10. A total of 500 eigenstates around the central one are 
considered; both a and /3 run from -250 to 250. 

In Fig.|9]we show the matrix elements for n{k = 0) and 
N{k = tt). We fix an effective temperature T = 5 and pick 
the eigenstate that has energy closest to it as the central state 
and 500 states more around it. These 501 eigenstates are used 
to compute the matrix elements, where a and /3 correspond to 
each one of the states, running from -250 to 250. Overall, the 
further away the element is from the diagonal, the smaller it 
becomes. The off-diagonal elements can be seen to be very 
small in the chaotic regime, so one expects the time fluctua- 
tions after relaxation to be small. At integrability, V' ~ 0, 
the off-diagonal matrix elements are seen to be slightly larger 
than in the chaotic regime, but it is only for large values of 
V' where we find that the off-diagonal elements become very 
large. So, in the latter regime, as expected, time fluctuations 
after relaxation will be lar ge. This is in agreement with the 
dynamics observed in Ref. 123]. 

The results for the matrix elements of K and / (Fig.fTOl) are 
qualitatively similar as the ones for h{k ~ 0) and N{k = tt) 
(Fig. |9]l. However, quantitative differences are also evident 
and hint that the dynamics of different observables in experi- 
ments may exhibit quantitative differences, particularly when 
the system is approaching localization for large values of V'. 
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FIG. 10: (Color online.) Matrix elements of K [top panels] and I 
[bottom panels]; L = 22, A: = 0. The data correspond to the same 
eigenstates of the Hamiltonian as in Fig.|9] In the bottom right panel 
(for V' — 10), the diagonal elements of / are beyond the interval 
presented in the plot. 

VI. CONCLUSIONS 

We have studied half-filled quantum chains of hard-core 
bosons with repulsive interactions. The considered systems 
are integrable in the presence of nearest-neighbor (NN) hop- 
ping and interactions. The addition of next-nearest-neighbor 
(NNN) interactions, which are characterized by the parame- 
ter V', leads to two different transitions: from integrability to 
chaos, as V increases from zero, and then from chaos to lo- 
calization in momentum space, as y — s- oo and the system 
approaches the atomic limit. We have investigated the valid- 
ity of the eigenstate thermalization hypothesis (ETH) in the 
three regimes. ETH is found to hold whenever chaos devel- 
ops. In this domain, the eigenstate expectation values (EEVs) 
for states close in energy become very similar. 

Our results have confirmed previous works ll22l l27ll which 
showed that quantum chaos in finite systems, and thus the va- 
lidity of ETH, depends on the system size L and the range of 
interactions. As L increases, the transition to chaos happens 
to smaller values of V', but the extrapolation to the thermo- 
dynamic limit still requires further studies. In terms of inter- 
actions, the Hamiltonian describing our system is a banded 
matrix, since it is restricted to two-body interactions. As a re- 
sult, chaos develops only in the center of the spectrum; at the 
edges, the eigenstates are more localized. Initial states with 
energy close to the borders of the spectrum may therefore be 
unable to thermalize. 

In the present work, we have focused on how the onset of 
chaos and the behavior of the EEVs may be affected by two 
other factors: the transition to localization and the presence 
of symmetries. On the way, we have shown that the opening 
of a gap in the ground state does not prevent thermalization. 
Even if the ground state of the system becomes an insulator. 



the structures of the eigenstates in the chaotic domain and 
close in energy, as well as the corresponding EEVs, do not 
fluctuate away from the edges of the spectrum. All EEVs of 
eigenstates close in energy were seen to become very similar, 
independently of the /j-sector they belong to and of discrete 
symmetries that may not have been accounted for during the 
diagonalization. This means that we do not find signs of rare 
states [29] in those systems. Moreover, the range of values of 
V' over which the ETH holds increases with L, carrying the 
viability of thermalization deeper into the insulating phase. 
This corroborates our findings in Ref. IztIi . where a system 
with 1/3 filling was considered. 

In addition to the conservation of the total number of parti- 
cles, the systems that we have analyzed presented also transla- 
tional, reflection, and particle-hole symmetries. Chaos indica- 
tors based on the eigenvalues may miss the transition to chaos 
when different symmetry sectors are mixed. Level repulsion 
is a main feature of chaotic systems, but to be noticeable it 
requires the examination of each symmetry sector separately. 
This is not the case for quantities that depend on the eigenvec- 
tors. In the chaotic domain, eigenstates from different sym- 
metry sectors were still found to have very similar structures. 
Therefore, the lack of fluctuations of quantities measuring the 
complexity of the eigenvectors, such as the structural entropy, 
or similarly, the lack of fluctuations of EEVs for eigenstates 
close in energy, has been shown to be a reliable alternative to 
identifying the chaotic region, especially in situations where 
unknown symmetries may be present. 

In terms of what to expect for the dynamics, we have shown 
that in the chaotic regime the off-diagonal elements of the 
few-body observables of interest are very small, so that time 
fluctuations after relaxation are expected to be small. At in- 
tegrability, V' = 0, off-diagonal elements were found to be 
slightly larger than in the chaotic regime, but it was only for 
large values of V', when the system approaches localization, 
that very large off-diagonal matrix elements were seen. As 
expected, in the latter regime the relaxation dynamics will be 
very slow and time fluctuations after relaxation will be large. 
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Appendix A: Eigenstates and Observables 

The purpose of this appendix is to provide further illustra- 
tions for the structure of the eigenvectors and for the EEVS 
across the two transitions achieved by increasing V', from in- 
tegrability to chaos and from chaos to locaUzation in the mo- 
mentum basis. 
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1. Shannon entropy 

Figure [TT] shows the Shannon entropy ^ in the fc-basis for 
various values of V' . The results are comparable to those for 
IPR in Fig. [T] becomes a smooth function of energy in 
the chaotic regime [for L = 22, when < y < 5]. Here, 
two curves are distinguished in the middle of the spectrum. 
The lower curve is associated with the more limited capabil- 
ities for spreading of the eigenstates in sectors k = and 
k = L/2, where particle-hole symmetry is also present. It 
is remarkable, that even when all fc-sectors are combined to- 
gether, the Shannon entropy is still capable of identifying the 
chaotic region. 

For V' 0, the eigenstates have very different levels of 
delocalization, even when close in energy. For V' oo, 
the eigenstates divide into energy bands and localize in the 
/c -basis (small values of Sfe). 
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FIG. 1 1 : (Color online.) Shannon entropy in the fc-basis vs energy 
per site for all fc-sectors; L — 22. Black dots: sectors fc = and 
k = L/2; light (red) dots: all other fc-sectors. 

Figure[T2]depicts similar results as Fig.[TT]but for a smaller 
system, with L = 20. The comparison between Figs. [TT] 
and [12] clearly shows that the fluctuations of S^ reduce in the 
chaotic regime and the separation between Sk for k ~ 0, -L/2 
and Sk for the other fc-sectors increases with L. This gap may 
also be taken as an indication of the chaotic regime. With in- 
creasing system size, the chaotic regime starts at smaller val- 
ues of V' and moves towards larger values of V' in the region 
where localization in fc-space starts to become evident by the 
reduction of the values of Sfc. 



2. Effect of AE in the Fluctuations of the EEVs 

Figure [T3] shows the same results from Fig. [6] for L = 22, 
but now for different values of A-E. The results are not much 
affected by the exact value of AE and the overall behavior is 
still the same: the EEVs approach the microcanonical average 



FIG. 12: (Color online.) Shannon entropy in the fc-basis vs energy 
for all fc-sectors; L = 20. Black dots: sectors fc = and k — L/2; 
light (red) dots: all other fc-sectors. 



in the chaotic region, but show large fluctuations in the inte- 
grable and localization regimes. We notice that in our plots, 
the larger relative fluctuations and larger effects of the win- 
dow of energy AE, which occur for the kinetic energy and 
the interaction energy in the chaotic regime, are related to the 
fact that their EEVs (see Fig.|4]i, and hence their mean values, 
approach zero for the windows of eigenstates selected for our 
calculations. 




FIG. 13: (Color online.) Average relative deviation of eigenstate 
expectation values with respect to the microcanonical result as a 
function of V' ; L = 22. The average is performed over all eigen- 
states (including all momentum sectors) that lie within the window 
[E — AE, E+AE]. The ten curves shown in each panel are obtained 
for AE — 0.05, 0.1, 0.15, . . . 0.5; the thick dashed line corresponds 
to AE = 0.4, which is the value considered in Fig.[6] The effective 
temperature is T = 5. 
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